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ABSTRACT 

We examine the possibility of the active galactic nucleus (AGN) obscuration by dusty gas clouds that 
spurt out from circumnuclear starburst regions. For the purpose, the dynamical evolution of gas clouds 
is pursued, including the effects of radiation forces by an AGN as well as a starburst. Here, we solve 
the radiative transfer equations for clouds, taking into consideration the growth of clouds by inelastic 
cloud-cloud collisions and the resultant change in optical depth. As a result, it is shown that if the 
starburst is more luminous than the AGN, gas clouds are distributed extensively above a galactic disk 
with the assistance of radiation pressure from the starburst. The total covering factor of gas clouds 
reaches a maximum of around 20%. After several lO^yr, gas clouds with larger optical depth form by 
cloud-cloud collisions and thereafter the clouds fall back due to weakened radiation pressure. The larger 
clouds undergo runaway growth and are eventually distributed around the equatorial plane on the inner 
sides of circumnuclear starburst regions. These clouds have an optical depth of several tens. The result 
is qualitatively consistent with the putative tendency that Seyfert 2 galaxies appear more frequently 
associated with starbursts than Seyfert Is. On the other hand, if the AGN luminosity overwhelms that 
of the starburst, almost all clouds are ejected from the galaxy due to the radiation pressure from the 
AGN, resulting in the formation of a quasar-like object. The origin of obscuration of AGNs is discussed 
with relevant observations. 
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1. INTRODUCTION 



The origin of the active galactic nucleus (AGN) ob- 
scuration is a key issue for a thorough understanding of 
AGNs. In the context of the unified model, an obscur- 
ing torus has been thought to be responsible for the AGN 
obscuration, and the dichotomy of AGN types has been 
attributed to the viewing angle toward the AGN (see An- 
tonucci 1993 for a review). However, recent observations 
on circumnuclear regions of AGNs have gradually revealed 
that Seyfert 2 galaxies appear to be more frequently as- 
sociated with starbursts than Seyfert Is (Heckman et al. 
1989; Maiohno et al. 1997, 1998a; Perez-Olea & Cohna 
1996; Hunt et al. 1997; Malkan, Gorjian, & Tam 1998; 
Schmitt, Storchi-Bergmann, & Fernandes 1999; Storchi- 
Bergmann et al. 2000; Gonzalez Delgado, Heckman, & 
Leithere 2001). Furthermore, quasars are mostly observed 
as type 1 AGNs, regardless of the star formation activity 
in host galaxies (Barvainis, Antonucci, & Coleman 1992; 
Ohta et al. 1996; Omont et al. 1996; Schinnerer, Eckart, & 
Tacconi 1998; Brotherton et al. 1999; Canalizo & Stockton 
2000a, 2000b; Dietrich & Wilhelm-Erkens 2000; Solomon 
et al. 2004). In addition, Ueda et al. (2003) have recently 
found from the hard X-ray luminosity function of AGNs 
that the fraction of X-ray absorbed AGNs decreases with 
the AGN luminosity. These recent findings seem beyond 
understanding in the context the unified model. 

The frequent association of circumnuclear starbursts 
with Seyfert 2 galaxies suggests that the obscuring of 
AGNs could be physically related to the starburst events. 
Ohsuga & Umemura (1999, 2001) have considered a phys- 
ical mechanism to link the AGN type to the circumnuclear 



starburst, where an obscuring wall composed of dusty gas 
can form with the assistance of the radiation pressure from 
circumnuclear starbursts. More recently, Wada & Norman 
(2002) have performed three-dimensional hydrodynamic 
simulations on the nuclear starburst regions, and found 
that the gas blown out from a galactic disk due to multiple 
supernova explosions in a starburst builds up a torus-like 
structure in a highly inhomogeneous and turbulent man- 
ner. Both models may provide potential mechanisms to 
relate the starburst events with the obscuring of AGNs. 
But, Ohsuga & Umemura focused just on the equilibrium 
configuration of dusty gas, and Wada & Norman did not 
include radiation forces, which are likely to play an impor- 
tant role in the early phase of starbursts. Hence, in order 
to construct a more realistic model, the dynamical evo- 
lution of obscuring materials should be solved, including 
radiation forces. This requires a radiation-hydrodynamic 
(RHD) simulation in three-dimensional space. However, a 
full RHD simulation is extremely time consuming even us- 
ing the state-of-the-art computational facilities. In this 
paper, we solve the dynamics of discrete clouds under 
the adhesion approximation, where the growth of clouds 
is brought by inelastic cloud-cloud collisions. Simultane- 
ously, to properly include the radiation forces, we solve the 
radiation transfer equation for clouds, taking the change 
in optical depth due to the cloud growth into account. In 
§2, the circumnuclear regions and gas clouds are modeled. 
In §3, we formulate the equation of motion for gas clouds 
including radiation forces. In §4, the numerical results are 
shown, and the distribution of gas clouds and the obscu- 
ration of AGN are analyzed. In §5, based on the present 
results, we discuss the origin of AGN obscuration and the 
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relation to narrow line regions. §6 is devoted to the con- 
clusions. 

2. MODEL 

Recent high-resolution observations have revealed cir- 
cumnuclear starburst regions with a radial extension of 10 
pc up to 1 kpc, which frequently exhibit a ring-like fea- 
ture (Pogge 1989; Wilson et al. 1991; Forbes et al. 1994; 
Marconi et al. 1994; Mauder et al. 1994; Buta, Purcell, 
& Crocker 1995; Barth et al. 1995; Leitherer et al. 1996; 
Maoz et al. 1996; Storchi-Bergmann, Wilson, & Baldmin 
1996; Elmouttie et al. 1998). In this paper, we suppose 
a starburst region in ring configuration with a radius of 
200 pc and the total mass of starburst ring is assumed 
to be MsB — 1O^M0, similar to the model by Ohsuga & 
Umemura (2001). Also, we settle a galactic bulge compo- 
nent with the mass of Mqb = lO^^M© and the radius of 
iZcB = 1 kpc. 

In some galaxies, a galactic bulge solely cannot explain 
the observed rotation curve. For instance, in the circum- 
nuclear regions of the Circinus galaxy, the rotation velocity 
requires an additional component with 10'^^ ^"M0 within 
a starburst ring (Elmouttie et al. 1998). Also, the stel- 
lar velocity dispersion suggests that the Circinus galaxy 
has a spread component within a radius of several 100 pc 
(Maiolino et al. 1998a). In this paper, with taking such 
mass distributions into consideration, we also assume an 
inner bulge component with the mass of Myb = IO^Mq 
and the radius of i?iB = 100 pc. As for a central su- 
permassive black hole (SMBH), the mass is assessed in 
terms of the recently inferred black hole-to-bulge mass re- 
lation, that is, A^BH A'/buigc ~ lO^'' (Richstone et al. 1998; 
Marconi & Hunt 2003, and references therein). Here, the 
SMBH mass is set to be Mbh = IO^Mq. 

We assume the energy spectrum of AGN to be in the 
form of ^AGN °^ between 0.01 and 100 keV (Bland- 
ford et al. 1990). The bolometric luminosity, Lagn, is 
set to be constant for a period of 10^ yr, which is a typ- 
ical age of AGN (basically Eddington timescale). In the 
starburst regions, we assume a Salpeter-type stellar initial 
mass function (IMF) for a mass range of [TOi,mu]; 

(/.(m,) =<^o(m*/MQ)-(i+^), (1) 
where s = 1.35 and is the stellar mass. Although 
mi and rriu in starburst regions are under debate, some 
authors argue that the IMF in a starburst is deficient in 
low-mass stars, where m\ « 2M0 is derived. The up- 
per limit is inferred to be « 4OM0 (Doyon, Puxley, 
& Joseph 1992; Chariot et al. 1993; Doane & Mathews 
1993; Hill et al. 1994; Blandle et al. 1996). We assume 
these values for m\ and rriu. The star formation rate, M*, 
is controlled by M* = (Msb/^sf) exp(— i/tsp), where t 
is the elapsed time after the initial starburst and tsF is 
the duration of star formation, which is assumed to be 
lO^'yr (e.g. Efstathiou, Rowan-Robinson, & Siebenmorgen 
2000). We employ a stellar mass-to-luminosity relation as 
{U/Lq) = [rriit/MQy with q = 3.7, and a stellar mass-to- 
age relation as U = 1-1 x lO^°yr(m*/M0)'^ with cj = 2.7 
(Lang 1974). All stars are assumed to emit the radiation 
in the blackbody. 

Next, we model the mass ejection from the starburst 
regions. The mass ejection is driven by multiple super- 
nova (SN) explosions in a starburst. Stars heavier than 



Wcrit = 8M0 are destined to undergo Type II SN explo- 
sions. The SNe restore almost all the mass into interstellar 
space and release the energy with a conversion efficiency to 
the rest mass energy as e « 10~^. Thus, the rate of mass 
restoration by SNe is assessed by AfsN(i) = isN(^)/c^£, 
where MsN(i) is the total mass of SNe that explode by 
time t and isN(i) is the total supernova luminosity that 
is temporally averaged. The energy input by multiple SN 
explosions leads to the formation of an elongated hot cav- 
ity (superbubble) , and consequently the shock- heated inter- 
stellar gas spurts out from the starburst regions (Shapiro & 
Field 1976; Tomisaka fc Ikeuchi 1986; Norman & Ikeuchi 
1989). A recent simulation by Wada & Norman (2002) 
have shown that the shock-heated gas eventually frag- 
ments into cold, dense gas clouds owing to the cooling. 
Hence, we simply suppose that gas clouds spurt from the 
starburst regions with the mass restoration rate estimated 
above. (It is noted that a part of preexisting interstel- 
lar gas could be also accelerated by the superbubble and 
therefore participate in the mass ejection from a starburst. 
Hence, the present mass ejection rate is regarded as a min- 
imum rate.) On the basis of Wada & Norman's results, we 
set the initial cloud radius and mass to be rc = 3.5 pc and 
Mc = IO^Mq as a fiducial model. Then, the hydrogen 
density in a cloud is 2 x 10^ cm~^. We also examine the 
effect of a different initial size of clouds by setting rc = 2.5 
pc or Tc — 4.5 pc. The cloud velocity is assumed to be 
Vc ^ 150 km/s, which is on the order of escape velocity. 
Figure 1 shows the evolution of mass ejection rate (solid 
line) and luminosity (dashed line) in the starburst regions. 
Before 4.3 x 10^ yr, the starburst is super-Eddington lu- 
minous for gas clouds of initial size. As for the effects 
of photoionization and photodissociation by the radiation 
from a starburst as well as an AGN, Ohsuga & Umemura 
(2001) have analyzed the inner structure of an optically- 
thick gas slab by solving the force balance and the energy 
equation. They found that the dust cooling is still effec- 
tive for a gas slab located at several hundred parsecs from 
the galactic center, and the resultant gas temperature be- 
comes on the order of 100 K. Hence, we assume the cloud 
temperature to be 100 K and the optical depths of clouds 
to be determined by the dust opacity. 

The spurted clouds can collide with each other in the 
space above the galactic disk. Here, we simply assume 
that an inelastic collision occurs if the distance between 
cloud centers is less than the sum of cloud radii. Also, it 
is assumed that the density of clouds does not alter after 
the collision, so that a cloud of larger size is generated by 
the collision. Taking into account the resultant change in 
optical depth of clouds, we calculate the dilution of the ra- 
diation from a starburst and an AGN. It should be noted 
that this assumption on the collision is quite naive. Naga- 
sawa & Miyama (1987) have simulated the head-on colli- 
sions of isothermal clouds, using a 3D hydrodynamic code, 
and have found that the cloud-cloud collisions are very in- 
elastic even for the Mach number of 10, and they result in 
the coalescence without disruption. However, in our sim- 
ulation, the Mach number can be up to several 10 or ^ 
100. Also, the off-center collisions can frequently occur. 
In this case, the cloud-cloud collision may result in not 
only the coalescence but also the disruption of clouds or 
the induction of star formation. These effects seem worth 
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incorporating in the future simulations. 

The self-gravity of clouds is ~ 10~^ of the external grav- 
ity in the model galaxy. Hence, the self-gravity of clouds 
can be neglected. But, the growth of clouds induced by 
collisions can lead to the gravitational instability, if the 
cloud radii exceed the Jeans scale, which is 10.8 pc for 
the assumed cloud temperature of 100 K. In the present 
simulation, we treat a cloud above the Jeans scale as a 
collapsing cloud which contributes no more to the obscu- 
ration. 

3. BASIC EQUATIONS 

The equation of motion is given by 



/rad.SB + /rad.AGN + /grav,SB + /grav,BH + /grav, bulge, 

(2) 

where r is the position vector of a cloud from the galactic 

center, and /grav.se, /grav,BH, and /grav, bulge are the gravi- 
tational forces from the starburst ring, the central SMBH, 
and the bulges, respectively. The radiation pressure force 
from the starburst, /rad,SB, is given by 



Jrad.SB = / ^——^ — -LgB 



4tt\1 



SB I 



1 - exp(-TSB) 
TSB 



(3) 



where xsb is the mass extinction coefficient, which is 
averaged over the spectrum of the starburst; xsb = 
J^XvLg^df/LsB, where is the mass extinction coef- 
ficient and is the spectrum. Since the mass den- 
sity is dominated by gas with little contribution of dust 
and the Thomson scattering is negligible in the opacity 
(Umemura, Fukue, & Mincshige 1997, 1998), Xi^ is given 
by = i^t/Pg' where k'^ is the opacity by dust grains 
per unit volume and pg is the gas density in a cloud, tsb 
is the cloud optical depth estimated by xsb, PSB{t) is the 
luminosity density in the starburst region, and ^sb is the 
position vector of a cloud from a volume element dV in 
the starburst region. We suppose the dust-to-gas mass 
ratio to be 0.03, which is three times as high as that ob- 
served in the solar neighborhood, because the metallic- 
ity in QSOs is found to be several to ten times higher 
than the solar metallicity (e.g. Hamann & Ferland 1993). 
Here, we employ a grain size distribution of a power-law as 
nd(ad) cc a^^'^ in the range of 0.01-1 /xm, which is found 
in the interstellar matter (Mathis, Rumpl, & Nordsieck 
1977), and the absorption cross section is determined by 
7ra^min[l, (27rz^ad/c)^], where is the grain radius. The 
density of solid material within a grain is assumed to be 
1.0 g cm~^ (e.g. Spitzer 1978). We find the resultant 
optical depth of a gas cloud to be Tc ~ 4. 

Similarly to the starburst case, the radiation pressure 
force from the AGN, /rad,AGN, is calculated by 



? Xagn ^agn ^ 

/rad,AGN = — 

c 47r \r\^^ 



1 - exp(-rAGN) 



TAGN 



(4) 



where xagn is the mass extinction coefficient, which is av- 
eraged over the spectrum of the AGN, tagn is the optical 
depth estimated by xagn- Here, we neglect the radiation 
force from SNe, because the average luminosity of SNe 
in the assumed IMF is less than 10% of the total stellar 
luminosity of star bursts. 



4. NUMERICAL RESULTS 

Here, we show the numerical results for a starburst- 
dominant case and an AGN-doininant c;ase. In each case, 
the time variations of the spatial distributions of gas clouds 
are presented, and the obscuration of AGN by the clouds 
is analyzed. 

4.1. Starburst-Dominant Case 

We consider a case with -Lagn = lO^^L©. Then the 
AGN luminosity is lower than the starburst luminosity un- 
til ~ 3 x lO'' yr. To show the cloud distributions, we use 
the coordinates depicted in Figure 2. The resultant cloud 
distributions projected onto the rz— and the a;y— plane 
are shown in the left panels in Figure 3. The colors of 
points stand for the cloud radii, where blue points are 3.5 
pc, which corresponds to the clouds missing collisions, and 
red ones are the maximum radius of 10.8 pc, which is the 
Jeans radius. Using the resultant cloud distributions, we 
calculate the optical depth, r, toward the direction with 
9 and (p. Taking the impact parameter into account, r is 
given by 



r = 



E 

i,b{i)<rc{i) 



2 1- 



Tc(«), 



(5) 



where i denotes each gas cloud, b{i) is the impact parame- 
ter (the distance from the cloud center to the linc-of-sight), 
Tc(i) is the optical depth of the cloud along the diameter, 
and rc(i) is the cloud radius. If b{i) < rdi), the clouds con- 
tribute to the optical depth. The right panels in Figure 3 
show T for the each direction {6, 0), using the Mollweide's 
projection for the all-sky plot. The colors represent the 
levels of T, ranging from (white) to 50 (red). 

As shown in Figure 3, gas clouds are distributed ex- 
tensively above the starburst ring at the stage of ~ 10^ 
yr. The radiation pressure from the starburst works to ac- 
celerate clouds, because the starburst is super-Eddington 
luminous until 4.3 x 10^ yr for gas clouds of initial size. 
But, clouds are not distributed around the polar regions. 
This is because clouds are ejected from the rotating star- 
burst ring and thus possess angular momenta. At this 
stage, the formation of larger clouds by cloud-cloud col- 
lisions is not promoted to a great extent. At ~ 2 x 10^ 
yr, clouds with large optical depths increase by repetitive 
cloud-cloud collisions and they fall back due to the re- 
duced radiation pressure per unit mass. The larger clouds 
undergo runaway growth and are eventually distributed 
around the equatorial plane on the inner sides of circum- 
nuclear starburst regions. These clouds have an optical 
depth of several tens. The angular momentum transfer by 
cloud collisions near the equatorial plane is not so efficient 
to allow the effective mass accretion on to the AGN. There- 
fore, the large clouds stay there and can contribute to the 
obscuration of the AGN. At ~ 10^ yr, the starburst lumi- 
nosity becomes sub-Eddington for all clouds, and therefore 
gas clouds fall back one after another. In this stage, a lot 
of clouds exceed the Jeans size, and thus the number of 
obscuring clouds decreases rapidly. In Figure 4, the tran- 
sition of the mass spectrum of the clouds is shown. This 
figure demonstrates the collision-induced growth of clouds 
and the decrease of obscuring clouds by the Jeans insta- 
bility. 
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Next, we evaluate the covering factor of obscuring 
clouds. Figure 5 shows the time variations of the cov- 
ering factor, which is defined by the fraction in the all sky 
of the area with the optical depth (t) larger than an op- 
tical depth given in the abscissa. If we see the covering 
factor of T > 2 — 3, it changes from ~ 13% at 10'' yr to 
a maximum of^ 21% at 2 x lO'' yr. This result shows 
that the starburst-origin dusty clouds can contribute to 
the obscuration of AGN to a notable extent. After that, 
the covering factor lessens clue to the fall back of clouds. 

To see the effect by a different choice of the initial size 
of gas clouds, we also calculate the models with the initial 
cloud radius of = 2.5 or 4.5 pc without changing the 
cloud mass. In Figure 6, we show the resultant covering 
factor against the optical depth. If the spatial distribution 
is the same, the shrink of clouds would result in the reduc- 
tion of covering factor due to the smaller geometrical cross 
section of clouds. The reduction of covering factor would 
be a factor of (2.5/3.5)^ = 0.51, if = 2.5 pc is assumed. 
But, as seen in Figure 6, the actual reduction of peak cov- 
ering factor is 0.16/0.21 = 0.76. This can be understood 
by the effect of radiation pressure. If the initial cloud size 
is smaller, the optical depth of a cloud becomes larger. 
Hence, the radiation pressure on a cloud is weakened, so 
that clouds are distributed in more compact space. This 
effect tends to enlarge the covering factor. On the other 
hand, for the case of Tc = 4.5 pc, the peak covering fac- 
tor increases by just a small factor of 0.23/0.21 = 1.1, 
although the geometrical cross section would lead to the 
increase of covering factor by (4.5/3.5)^ = 1.65. There is 
again the radiation pressure effect. Expanded clouds can 
be distributed more sparsely by less diluted radiation pres- 
sure, so that the covering factor is reduced. In addition, 
larger clouds more frequently collide with each other. This 
effect also reduces the covering factor. As a result, it turns 
out that the peak covering factor is not so sensitive to the 
initial size of clouds. 

4.2. AGN-Dominant Case 

Here, we show the results in an AGN-dominant case, 
where Lagn = 3 x IO^^Lq is assumed. In this case, the 
AGN luminosity is higher than the starburst luminosity 
all the time. Figures 7 and 8 show the results by the same 
quantities as Figures 3 and 5, respectively. The assumed 
AGN luminosity is close to the Eddington luminosity for 
the total mass of the galaxy. Hence, gas clouds are vig- 
orously accelerated. Even at the early stage of ~ lO'' yr, 
clouds are already distributed over kpc, as shown in Figure 
7. Very few collisions occur, because the radiation pres- 
sure from the AGN drives a spreading outflow of clouds. 
At ~ 2 X 10^ yr, a large fraction of clouds escape from 
the galaxy, so that most of gas clouds cannot contribute 
to the obscuration of AGN. At ~ 10^ yr, a small number 
of collision-induced large clouds are left. 

Due to the cloud motion driven by the AGN luminos- 
ity, the covering factor of obscuring clouds cannot reach a 
higher level. The maximum covering factor is ~ 5% at lO'^ 
yr, as shown in Figure 8. Hence, the AGN is likely to be 
observed as a type 1 with high probability. Such an AGN- 
dominant case may correspond to quasar events. The rela- 
tion between AGN type and luminosity is discussed below 
in the context of the obscuration. 



5. DISCUSSION 

5.1. Origin of Obscuration 

In the previous section, it is shown that the covering 
factor by starburst-origin dusty clouds increases in a lu- 
minous phase of the circumnuclear starburst. This can 
provide a qualitative explanation for the fact that Seyfert 
2 galaxies are more frequently associated with starbursts 
than Seyfert Is. On the other hand, obscuring clouds are 
not distributed around polar regions, as shown in Figure 
3. Thus, in the face-on view, the AGN is likely to be iden- 
tified as a type 1 even in a luminous phase of starburst, 
if there is no obscuring materials other than starburst- 
origin dusty clouds. In that sense, Seyfert Is may be 
coupled with starburst events. Interestingly, hidden star- 
bursts have been found in Seyfert 1 galaxies (Imanishi & 
Dudley 2000; Rodriguez- Ardila & Viegas 2003), and also 
there is an example of Seyfert 1 galaxies accompanied by 
circumnuclear starbursts, that is, NGC 7469 (Genzel et al. 
1995). In the present simulation, if we see the starburst 
ring at face-on view, starburst-origin dusty clouds obscure 
30% of the starburst ring at maximum. This covering fac- 
tor is not sufficient to produce perfectly hidden starbursts. 
However, further effects can change the covering factor. In 
the present simulation, we have ignored collapsing clouds 
above the Jeans scale. But, the collapsing clouds are likely 
to result in the star formation. Hence, the re-ejection of 
matter by supernovae in the collapsed object may also 
contribute to the obscuring. Furthermore, in the light of 
the tight correlation between SMBHs and galactic bulges 
(Richstone et al. 1998; Marconi & Hunt 2003; Kawakatu 
& Umemura 2004), a starburst in a bulge may contribute 
to the obscuration of nuclear regions (e.g. Umemura 2001; 
Kawakatu, Umemura, & Mori 2003). 

As for the ratio of Seyfert 2s to Is, Maiolino & Rieke 
(1995) concluded that Seyfert 2s appear to be four times 
more numerous than Seyfert Is. This implies that the to- 
tal covering factor of obscuring materials is likely to be 
of order 80%. In the present simulation, the maximum 
covering factor is around 20%. A hydrodynamic simula- 
tion by Wada & Norman (2002) also shows the covering 
factor to be 40% at an early phase (1.6Myr) of starburst. 
Hence, it seems that additional obscuring materials are 
required to account for the number ratio between Seyfert 
Is and 2s. Recently, we have some significant pieces of 
observational information about the obscuring materials. 
The Ay of the nuclear or circumnuclear regions is esti- 
mated to be between a few and several magnitudes by IR 
and optical observations (Rix et al. 1990; Roche et al. 
1991; Goodrich, Veilleux, & Hill 1994; McLeod & Rieke 
1995; Oliva, Marconi, & Moorwood 1999). On the other 
hand. X-ray observations have shown that most Seyfert 
2 nuclei are heavily obscured along the line of sight with 
at least Ay > 10 mag and sometimes Ay > 100 mag 
(Matt et al. 1996, 1999; Maiohno et al. 1998b; Bassani 
et al. 1999; Risaliti, Maiohno, & Salvati 1999). Also, it is 
argued that a component of obscuring materials must be 
extended up to > 100 pc in addition to a compact com- 
ponent confined to subparsec scales (Rudy, Cohen, & Ake 
1988; Miller, Goodrich, & Mathews 199l'; Goodrich 1995; 
McLeod & Rieke 1995; Maiolino et al. 1995; Maiolino & 
Rieke 1995; Malkan et al. 1998). These findings suggest 
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that the distributions of dusty gas around an AGN are 
much more diverse than previously considered. The dusty 
clouds ejected from starbursts may be the origin of the 
extended obscuring matter. In addition, inner obscuring 
materials (possibly a dusty torus) may cooperatively work 
to obscure the nucleus. The combination of outer and in- 
ner obscuring materials will be an important issue to be 
scrutinized. 

For the case that the AGN luminosity is predominant, 
dusty clouds arc blown out from a galaxy by the AGN ra- 
diation pressure, and therefore cannot contribute much to 
the obscuration. This result matches the fact that quasars 
are mostly observed as type 1 AGNs, regardless of the 
star formation activity in host galaxies (Barvainis, An- 
tonucci, & Coleman 1992; Ohta ct al. 1996; Omont ct al. 
1996; Schinnerer, Eckart, & Tacconi 1998; Brotherton et 
al. 1999; Canalizo & Stockton 2000a, 2000b; Dietrich & 
Wilhelm-Erkcns 2000; Solomon et al. 2004). Also, based 
on the hard X-ray luminosity function of AGNs, Ueda et 
al. (2003) have found that the fraction of X-ray absorbed 
AGNs decreases with the AGN luminosity. This trend is 
also consistent with the present picture of the radiation 
pressure-induced outflow. 

5.2. Relation to Narrow Line Regions 

It may be intriguing to consider the possible relation 

between dusty clouds and narrow line regions (NLRs) in 
AGNs, because the velocity dispersion of starburst-origin 
clouds, a few 100 km s~^ , is just on the order of that found 
in NLRs. 

One of common properties of NLRs is a fairly constant 
ionization parameter, i.e., U = S/nc ^ 0.01, where S is 
the ionizing photon flux and n is the density at irradi- 
ated cloud surfaces. Recently, Dopita et al. (2002) ar- 
gued that the ionization parameter of NLRs is success- 
fully accounted for by radiation pressure-dominated dusty 
photoablating clouds. This picture seems quite consistent 
with the present model. In addition, it is reported that 
Seyfert 1 galaxies have a smaller number of narrow-line 
clouds than Seyfert 2s (Schmitt 1998), and also the ex- 
tension of NLRs in AGNs is roughly proportional to the 
square root of the [O III] luminosity (Bennert et al. 2002). 
These trends are qualitatively compatible with the present 
results that the number of dusty clouds decreases with a 
dimming starburst, and the spatial extent of dusty clouds 
augments with increasing AGN luminosity. 

6. CONCLUSIONS 



We have investigated the dynamics of gas clouds that 
spurt out from circumnuclear starburst regions, including 
radiation forces from a starburst as well as an AGN. The 
results are summarized as follows: 

(1) In a case where the starburst luminosity is higher 
than the AGN luminosity, gas clouds are distributed at 
several hundred pc above the galactic disk, due to the 
radiation pressure from the starburst. The covering fac- 
tor of these clouds reaches around 20% at the maximum. 
According as the starbiu'st dims, gas clouds with large op- 
tical depth are formed by cloud-cloud collisions and fall 
back due to the reduced radiation pressure by starburst. 
The larger gas clouds which have optical depth of sev- 
eral tens are eventually distributed near to the equatorial 
plane. This starburst-dominant case is qualitatively con- 
sistent with the fact that Seyfert 2 galaxies appear to be 
more frequently associated with starbursts than Seyfert 
Is. 

(2) In a case where the AGN luminosity is dominant and 
super-Eddington, most gas clouds are blown out by the 
radiation pressure from the AGN. Resultantly, starburst- 
origin gas clouds can hardly contribute to the obscuration 
of the AGN. Hence, the AGN is likely to be observed as 
type 1 with high probability. This AGN-dominant case 
corresponds to a quasar-like object. Furthermore, it is 
predicted that the fraction of obscured AGNs decreases 
with increasing AGN luminosity. This prediction matches 
the fraction of X-ray absorbed AGNs found by Ueda et al. 
(2003). 

(3) The velocity dispersion of gas clouds is several 100 

km/s, which is just comparable to that found in NLRs. 
The observations for the ionization parameter of NLRs 
and the extension of NLRs depending on AGN luminosity 
are consistent with the present picture. 

In the present analysis, we have found that starburst- 
origin gas clouds can contribute to the obscuration espe- 
cially in large scales. Also, it has turned out that there 
is still something missing to account for all the properties 
of the AGN obscuration. It may be significant to consider 
further effects such as the combination of inner and outer 
obscuring materials. 
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Fig. 1. — The evolution of mass ejection rate {solid line) and luminosity of starburst region (dashed line). The starburst is super-Eddington 
luminous for gas clouds of initial size before the epoch indicated by a dot-dashed line (~ 4.3 X 10^ yr). 



z 




Fig. 2. — Schematic view of the circumnuclear structure and the coordinate system. 



http://www.rccp.tsukuba.ac.jp/Astro/watabe/ 

Fig. 3. — The distributions of gas clouds and optical depth to AGN. Lagn = lO^^L© is assumed. We show snapshots at 10^, 2 X lO'^, and 
10* yr. The left panels display the cloud distributions projected on the xy— and the r^— plane in the coordinates shown in Figure 2. The 
colors of points correspond to the radii of gas clouds. The right panels show t for the each direction (9, (p) using MoUweide's projection for 
the all-sky plot. The colors of points correspond to the levels of r. 
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Mo [10' M3] 

Fig. 4. — The mass spectrum of clouds at (a) 10^ yr, (b) 2 x 10^ yr, and (c) 10® yr. Horizontal axis is the cloud mass in unit of IO^Mq, 
while vertical axis is the cloud number. 
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Fig. 5. — The covering factor of the area whose line-of-sight optical depth is greater than an optical depth given in the horizontal axis. 
Lagn = is assumed. The evolution of the covering factor is shown at 1, 2, 3, 5, and 10 xlO'^ yr. 
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Fig. 6. — Same as Figure 5, but the initial cloud radius is 2.5 pc (above) and 4.5 pc (below) 
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Fig. 7. — Same as Figure 3, but Lagn = 3 x IO^^Lq is assumed. In this case, since the radiation pressure from AGN is high, most of gas 
clouds cannot contribute to the obscuration of AGN. 
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Fig. 8. — Same as Figure 5, but Lagn 



= 3 X 1O"L0. 



